Crystallization and gelation in colloidal systems with short-ranged attractive 

interactions 
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We systematically study the relationship between equilibrium and non-equilibrium phase dia- 
grams of a system of short-ranged attractive colloids. Using Monte Carlo and Brownian dynamics 
simulations we find a window of enhanced crystallization that is limited at high interaction strength 
by a slowing down of the dynamics and at low interaction strength by the high nucleation barrier. 
We find that the crystallization is enhanced by the metastable gas-liquid binodal by means of a 
two-stage crystallization process. First, the formation of a dense liquid is observed and second the 
crystal nucleates within the dense fluid. In addition, we find at low colloid packing fractions a fluid of 
clusters, and at higher colloid packing fractions a percolating network due to an arrested gas-liquid 
phase separation that we identify with gelation. We find that this arrest is due to crystallization 
at low interaction energy and it is caused by a slowing down of the dynamics at high interaction 
strength. Likewise, we observe that the clusters which are formed at low colloid packing fractions 
are crystalline at low interaction energy, but glassy at high interaction energy. The clusters coalesce 
upon encounter. 
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I. INTRODUCTION 

Crystallization, vitrification and gelation in systems 
with short-ranged attractive interactions like colloidal 
particles and proteins are widely studied phenomena P, 
0) S 0) IE B S B @- Crystallization kinetics is es- 
pecially important for proteins because good crystals 
are necessary for the determination of their 3D struc- 
ture. Furthermore, crystallization of biologically ac- 
tive proteins can lead to diseases like the human ge- 
netic catara ct fipj - While the crystal is an equilib- 
rium phase [ll| BQH, gels and glasses are 
non-equilibrium pha ses that occur due to dynamical ar- 
rest [S El H M S3- Vitrification in a fluid of N c 
hard spheres with diameter a c in a volume V , occurs for 



packing fractions rj — (Tra^/6)N C /V > 0.58 [16}, while 
the equilibrium phase diagram predicts a stable crys- 
tal phase for rj > 0.545 [2l|. In the repulsive colloidal 
glass phase, the particles are arrested as they become lo- 
calized in a "cage" formed by the surrounding spheres. 
Interestingly, the addition of non-adsorbing polymer, re- 
sulting in an effective attraction between the colloids, 
can enhance the relaxation process so that the system 
can crystallize. Nevertheless, for increasing attraction 
strengths the particles become more and more localized 
due to long-lived bonds between the neighboring parti- 
cles until the system is arrested in an attractive glass 
at smaller packing fractions [l9| . An arrested state at 



low rj is often called a gel. However, the relation be- 
tween gels and attractive glasses is still highly debated. 
In general, gels can be classified as chemical, if the par- 
ticles are linked by irreversible bonds, or as physical, if 
the bonds are not permanent. Physical gels can form 
through a non-equilibrium mechanism due to the arrest 
of a (metastable) gas-liquid separation by the attractive 
glass transition [22, [Hj]. On the other hand, equilibrium 
gelation can occur in systems where the gas-liquid phase 
separation is suppressed by e.g. long-ranged repulsions 
or directional interactions [8j . For increasing contact en- 
ergy these two mechanisms for the formation of physical 
gels continuously evolve towards a diffusion limited clus- 
ter aggregation (DLCA) that leads to chemical gels. The 
metastable gas-liquid separation is therefore important 
for the formation of non-equilibrium gels [22 |. and it is 
widely recognized as an important factor in protein crys- 
tallization [24], HE HE H3] • There has been much recent 
experimental, theoretical and simulation effort in under- 
standing the equilibrium phase diagram or the kinetics 
and dynamics of these systems, but only a few studies 
explored the relationship between the equilibrium phase 
diagram and the underlying kinetic pathways for gela tion 
and crystallization. Compared to previous works [28L |29| 
we systematically determine the crystallization behavior 
in relation to the equilibrium phase diagram in the re- 
gion of the metastable gas-liquid phase separation and 
we investigate the relation between crystallization and 
gelation. 
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II. MODEL AND METHOD 

We study a model system of N c hard spheres of diame- 
ter a c and non-adsorbing polymers with diameter a p de- 
scribed by the Asakura-Oosawa-Vrij (AOV) model [30l |. 
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The polymers induce an effective short-ran ged attractive 
depletion interaction between the colloids [30l l3lj . The 
total interaction potential is defined as 
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where r; 



\Ri—Rj I is the distance between two colloidal 



particles with Ri the center-of-mass of colloid i, and 
(3 = 1/ksT with ks the Boltzmann's constant and T the 
temperature. Thus the size ratio q — o p /a c determines 
the range of the interaction. The polymer reservoir pack- 
ing fraction r\ r v = {■na p /Q)p r p , is proportional to the poly- 
mer density p p in the reservoir with which our system is in 
osmotic contact. This parameter controls the strength of 
attraction. As a measure of this interaction strength we 
also use the contact energy (3U — ~(3Udcp{fij = o c ) and 
the experimentally accessible reduced second virial coef- 
ficient B* 2 = 1 + 3/o 3 c J^ +ap r 2 (l - exp(-(3U dev (r)))dr, 
which is normalized by the second virial coefficient of 
hard spheres with a diameter a^. Although colloids often 
have a Coulombic repulsion [32| and the attractive po- 
tential in proteins is often patchy [33j our simple model 
captures the basic phenomena that play a role in the com- 
petition between equilibrium and non-equilibrium phases 
of systems with short-ranged attractions. 

We determine the equilibrium phase diagram by calcu- 
lating the dimensionless free energy density / = (3F/V as 
a function of the colloid packing fraction r/ c and contact 
energy PU (polymer reservoir packing fraction r] p ) with 
Monte Carlo simulations. We use thermodynamic inte- 
gration to relate the free energy of the effective system 
to that of a reference system at the same colloid volume 
fraction rj c . To this end, we write the total free energy 
density as the sum of two contributions 



f(N c , V, r, r p ) = fc(N c , V) + f dcp (N c , V, r, r p ) 



(3) 



where / C (A C , V) is the free energy density of a system of 
N c hard spheres in a volume V, for which we used the 
Carnahan-Starling expression [34| for the fluid phase, and 
the Hall expression [35| for the face-centred-cubic (f.c.c.) 
solid phase. The free energy density /d cp is the contri- 
bution of the depletion potential @ to the free en ergy 
density, and it is computed using the A- integration [36j |. 
In order to map out the phase diagram we determine 
the total free energy density f(rj c , t£) for many state 
points (j] c ,r]p) in simulations, and we employ common 
tangent constructions at fixed jf p to obtain the coex- 
isting phases [13] • To trace the gas-solid binodal we 
use A c =512 particles with 15 x 10 4 Monte Carlo moves 
per particle for equilibration and 15 x 10 4 Monte Carlo 



moves per particle for production. To check the accu- 
racy of the results few simulations were repeated using 
15 x 10 5 Monte Carlo moves per particle for equilibration 
and 15 x 10 5 Monte Carlo moves per particle for pro- 
duction. To trace the metastable binodal while avoiding 
crystallization [53, [H, H§| we use a smaller system with 
N c =120 particles. To get good statistics 5 to 10 inde- 
pendent Monte Carlo runs are necessary for each value 
of the density. 

To study the dynamics of the system we perform 
standard Brownian dynamics (BD) simulations [4(| by 
quenching homogeneous configurations to one of the state 
points indicated by the points in Fig. [TJ The time step 
St is chosen to be St = 2 x 10~ 5 tb with tb — o"^/{ADq), 
where Do is the Stokes-Einstein diffusion coefficient. We 
neglect all hydrodynamic interactions between the parti- 
cles. The random forces mimic the interaction between 
particles and solvent, and are sampled from a Gaussian 
distribution with variance 2Do5t. The random force and 
the dissipative term provide the system with a heat bath 
at constant temperature. We carry out BD simulations 
using particle numbers N c ranging from 600 to 1200 and 
to avoid discontinuities in the potential we replace the 
hard-core potential with the equivalent repulsive soft po- 
tential 1/r 36 . We studied the evolution of a system for 
a total time t = 600tb ■ In order to study the dynamical 
arrest, we compute the bond correlation function 
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where riy=l if the separation between particle i and j is 
smaller than a cutoff value of 1.25 o c , otherwise ny=0. 
The cutoff value corresponds to the minima between the 
first and second peak of the radial distribution function. 
The function $s(t) is similar in spirit to the one used 
in [32j l. and counts which fraction of the bonds at t=0 is 
still present at time t. 

We characterise the local structure around a particle i 
by a set of numbers 
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where Yi m (fij) are spherical harmonics, fij is a unit vec- 
tor in the direction of the bond between particle i and 
particle j. The sum runs over all Nj,(i) neighbours of 
particle i. We then construct the dot product 



q(ij) = qi(i) ■ = E 1irn(i)q* m U) 



(6) 



m=—L 



where i and j are neighbouring particles, and q* m (i) is the 
complex conjugate of qi m (i). We normalized the vector 
qi(i), such that qi(i) ■ q}(i) — 1. In our analysis particles 
i and j are defined as neighbors if n,y = 1, i.e. Rij < 
1.25fi c . Two particles are defined to be joined by a crystal 
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bond if q(ij) > 0.5. Particle i is crystal-like if at least 7 of 
its bonds with neighbouring particles are crystal- like [4l| . 
As in other studies [HJ we use the 1=6 order parameter. 



III. RESULTS 

We determine the phase diagram by calculating free 
energies in Monte Carlo simulations [15j for the potential 
([1]) and for a size ratio q = 0.15, which corresponds ap- 
proximately to the experimental value of recent work [43| . 
Fig. [1] shows the equilibrium phase diagram in the con- 
tact energy f3U (polymer reservoir packing fraction r\ r v ), 
colloid packing fraction r\ c representation. For f3U =0 we 
recover the hard-sphere phase diagram with a fluid phase 
at rj c ~ 0.494 in coexistence with a face-centered-cubic 
(fee) crystal phase at rj c ~ 0.545. For increasing (3U the 
region of fluid-solid coexistence opens up. At the critical 
energy (3U cr ~3.48 ( {rfp) C r — 0.316) we find a gas-liquid 
transition which is metastable with respect to a broad 
fluid-solid coexistence. 

In order to follow the crystallization, we characterize 
the local structure around particle i by the q$ order pa- 
rameter [4l|. We define the crystal fraction f cr (t) as the 
number of crystal-like particles at time t divided by the 
total number of particles. The red (dark grey) region 
in Fig. [1] indicates the statepoints for which the crys- 
tal fraction is larger than 0.4 at t = 600tb while the 
orange (light grey) region indicates a fraction of crystal 
smaller than 0.4. Let us discuss first the low (3U limit 
of the red (dark grey) region. Clearly at low density 
the metastable gas-liquid binodal has a dramatic effect 
in inducing (spontaneous) crystallization of the colloidal 
particles, with the bottom boundary of the crystalliza- 
tion region as denoted by the red (dark grey) region in 
Fig. [1] coinciding with the metastable gas-liquid binodal. 
Below the binodal line, and at low colloid density, we do 
not observe spontaneous nucleation, in agreement with 
the observations of 29]. Neither do we observe, in the 
time scale of our simulations, crystallization just below to 
the critical point despite the presence of critical fluctua- 
tions that lower the nucleation barrier [44| • On the other 
hand, at medium-high colloid density we observe crys- 
tallization also outside the binodal. We now focus on 
the high f3U limit for crystallization. Fig. [5] shows that 
the crystal fraction slowly decays to zero for increasing 
attractive energy (3U, and that the crystal fraction does 
not depend strongly on rj c for fixed f3U. This indicates 
a crystallization process hindered by a slowing down of 
the dynamics. The boundary between the orange (light 
grey) and red (dark grey) regions in Fig. [1] is a horizon- 
tal line at (3U « 4.4. As it has been demonstrated for 
short-ranged attractive systems, isodiffusivity lines are 
also horizontal when approaching the gas-liquid binodal 
[1, HI]. Therefore, the slow decrease of crystallinity can 
be seen as an indicator of the proximity of the (attrac- 
tive) glass transition line that intersects the gas-liquid 
binodal. It is worthwile noting that the attractive glass 



transition was determined to be at ~ 0.43 using the 
long-time scaling of the incoherent correlation functions 
obtained from computer simulations for a system with 
a size ratio of q=0.1, and using a pair potential that in- 
cluded a long-ran ge re pulsive barrier to prevent gas-liquid 
phase separation [45|]. This value is remarkable close to 
our value of eta r p ~ 0.4 for which we find hardly any 
crystallization. 
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FIG. 1: (Color online) The equilibrium phase diagram of 
a model colloid-polymer mixture with size ratio q = 0.15 as 
a function of the colloid packing fraction rj c and interaction 
strength U/ksT (polymer reservoir packing fraction rf v and 
the reduced second virial coefficient B^)- The circles and 
the triangles indicate the stable fluid-solid coexistence region. 
(The circles are the results of simulations where a total of 
3 x 10 4 MC moves per particles were used. The triangles are 
the results of more accurate simulations where a total of 3 x 
10 5 MC moves per particles were used). The continuous lines 
are a guide to the eye. The squares indicate the calculated 
coexistence points of the gas-liquid binodal. The thin dashed 
line is a fit to the binodal. The points indicate the state points 
where we carried out BD simulations. The red (dark grey) 
region indicates where the crystal fraction is larger than 0.4 
at t = 600tb- The orange (light grey) region indicates where 
the crystal fraction is smaller than 0.4 at t = 600rs. 

Fig. [3] shows the phase diagram as in Fig.[T]but we la- 
bel state points that will be discussed below. To demon- 
strate the existence of a dynamical arrest we compute 
the bond correlation function <&B(t) (Fig. [J) measured 
starting at a time t=600 tb after the initial quench in 
the state points a-d in Fig. [3] We observe a divergence in 
the characteristic decay time between f3U = 3.3 (point b 
in Fig. [3J and (3U = 3.74 (point d in Fig. [3]), giving ev- 
idence of a fluid-solid (either crystalline or amorphous) 
transition. The behavior of ^^(i) is very similar whether 
the solidification is due to crystallization {(3U = 3.74) or 
an attractive glass transition ((3U = 6.6). 

We now focus our attention to the state points in- 
side the gas-liquid region. We find that in certain re- 
gions of the phase diagram the phase separation is ar- 
rested in a permanently percolated structure. The thick 
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FIG. 2: Crystal fraction /cr at t = 600rg as a function of the 
interaction energy U/ksT for (a) ?7 C =0.10. (b) r; c =0.25. (c) 
?? c =0.40. 




FIG. 3: (Color online) The same as Fig. [TJ but we indicate 
the state points discussed in the text. For clarity the binodal 
symbols are removed . The thick (blue) dashed line indicates 
the transition between percolating and non-percolating states. 

dashed line in Fig. [3J corresponds to the percolation 
limit. For quenches inside the gas- liquid coexistence we 
find, on the left side of the percolation line, clusters 
that are crystalline at low f3U (red (dashed) region in 
Fig. [3J and amorphous at high (3U (orange (light grey) 
region in Fig. [3j, while on the right-hand side, we find 
a 'long-lived' percolating network (gel) with high-density 
branches that are crystalline at low (3U (red (dashed) re- 
gion in Fig. [3]) and amorphous at high j3U (orange (light 
grey) region in Fig.[3j). 

We now analyze in more detail the kinetic pathways 
in different regions of the phase diagram. The first 
state point we consider is point e in Fig. [3] ( r? c =0.05, 
(3U = 3.74, and 7^=0.34), just beyond the gas binodal. 
The direct path to the equilibrium phase would be nu- 
cleation of the crystal phase in the dilute fluid phase. 
The actual kinetic pathway is shown in the attached 
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FIG. 4: (Color online) Bond correlation functions $s(t) as 
a function of time t for r] c —0.25. Circles indicate results for 
energy f3U — (point a in Fig.|3}, triangles-up indicate (3U = 
3.3 (point b in Fig. [3j, triangles-down indicate f3U = 3.74 
(point c in Fig. [3J), triangle-left indicate j3U = 6.6 (point d in 
Fig.0. 



movies |46j . Fig. EJa) and (b) show typical configurations 
at t = 50tb and t = 600tb, respectively. We observe first 
in Fig. [5fa) the formation of metastable liquid droplets 
(yellow/small light spheres) surrounded by gas particles 
(large light spheres). Subsequently the nucleation of the 
crystal (red/large dark spheres) starts within the fluid 
droplets. The crystalline clusters grow until a local equi- 
librium between the cluster and a surrounding gas phase 
is reached. These metastable clusters survive for a long 
time (Fig.[5Jb)) due to low diffusion, but irreversibly coa- 
lesce upon touching, in contrast with the findings of [43| . 
At higher colloid density r\ c = 0.25 and /3f7=3.74 (point 
c in Fig. [3]), we observe spinodal decomposition, and 
subsequently nucleation of crystal nuclei within the fluid 
branches of the spinodal structure. The crystalline clus- 
ters grow until the whole spinodal structure is crystalline. 
Fig.jnja) shows that the spinodal structure is arrested by 
crystallization. At higher interaction strength /3Z7 =6.6 
and ?7 C =0.25 (point d in Fig. [3J, we observe spinodal de- 
composition, and only very small crystal nuclei appear 
within the fluid branches of the spinodal structure. The 
spinodal decomposition is arrested into a branched gel 
by the formation of a glassy solid (Fig. (Hb)). This sce- 
nario is consistent with non-equilibrium gelation induced 
by spinodal decomposition arrested by an attractive glass 
transition. The metastable gas-liquid phase separation is 
arrested due to decreased diffusion of single particles by 
crystallization at low j3U. On the other hand, at high /3U 
the phase separation is arrested by a slow dynamics that 
is consistent with vitrification, but we do not explicitly 
calculate the exact location of the glass transition here. 

Fig. [7] shows the crystal fraction fci(t) as a function 
of time t for (3U—3.3, 3.52, 4.18, and 6.6 and for vary- 
ing rj c (filled symbols along the red lines A-D in Fig. [3J. 




FIG. 5: (Color online) Typical snapshots for the state point 
e in Fig. [3] (j7 c =0.05 and f3U = 3.74) at t=50 r s (a), and 
at t=600 tb (b). Gas particles are indicated by large light 
grey spheres, fluid particles are indicated by yellow (small 
light grey) spheres and crystal particles are indicated by blue 
(large dark grey) spheres. The complete kinetic pathway is 
shown in the attached movies [46l ]. 



The evolution of the crystal fraction is characterized by 
three distinct regimes. In the first regime the fraction of 
crystal is equal to zero. This regime corresponds to the 
time elapsed between the quench and the formation of 
a critical nucleus (induction time). The second regime 
corresponds to the fast growth of the crystal. Fig. [7fa) 
shows that for statepoints ?7 C =0.05, 0.20, and 0.40 of line 
A in Fig. [3J the induction time is significantly different 
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FIG. 6: (Color online) A snapshot for state point c in Fig. [3] 
(t? c =0.25 and fill = 3.74) at t=600 tb) (a) and for state 
point d in Fig. (^=0.25 and f3U = 6.6) at t=600 r B (b). 
Gas particles are indicated by large light grey spheres, fluid 
particles are indicated by yellow (small light grey) spheres 
and crystal particles are indicated by blue (large dark grey) 
spheres. The complete kinetic pathway is shown in the at- 
tached movies [46l ]. 



from zero. We note here that the induction time is larger 
than 600 tb for r\ c = 0.05 and 0.20. All these statepoints 
are outside the metastable gas-liquid binodal region. In 
Fig. [ZUb), we find that statepoint r\ c — 0.05 of line B 
in Fig. [3] shows no crystallization at all, and hence the 
induction time is again larger than 600t£. This state- 
point also lies outside the gas-liquid coexistence region, 
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FIG. 7: (Color online) Crystal fraction fci as a function of 
time t for (a) ,3(7=3.3 (line A in Fig.[3), (b) ,3(7=3.52 (line in B 
Fig. El), (c) £17=4.18 (line in C in Fig-EJ, (d) ,8(7=6.6 (line D 
Fig. [3}. Squares indicate a colloid packing fraction ^=0.05, 
circles indicate 77 c =0.20, triangles-up indicate r/ c =QA0, and 
triangles-down indicate rj c =0A5 (same symbols as in Fig. 3) 



i.e, rj c < 7y* ot , where r)l at denotes the coexisting packing 
fraction of the metastable gas phase. For quenches to 
r) c > rf c at , i.e., ry c =0.20, 0.40 and 0.45, we find that the 
induction time is close to zero. Similarly, for quenches 
inside the gas-liquid spinodal as shown in Fig. [T^c) and 
corresponding to line C and D in Fig. O the induc- 
tion time is also close to zero as the formation of the 
spinodal structure is followed immediately by nucleation 
of small crystallites and subsequent growth of the crys- 
tals. In this case, the induction time is related to the 
spinodal decomposition ~ td [H| that is not visible on 
the timescale of figure [7] To summarize, we find that the 
induction time is close to zero for all quenches beyond 
the gas binodal of the metastable gas-liquid coexistence 
region, i.e., for r\ c > rf^ 1 . The fraction of crystal grows 
rapidly after the induction time due to the huge differ- 
ence in chemical potential between the liquid and the 
growing crystal phase. In the third regime, the crystal 
fraction grows slowly by coarsening and defect removal 
processes. The fraction of crystal reached by the system 
depends strongly on the contact energy. At low f3U par- 
ticle diffusion is fast enough for crystallization. On the 
other hand, at high j3U the particles diffuse slowly and 
despite the thermodynamic drive to form the crystal the 
system is trapped in an amorphous solid. To summarize, 
for state points inside the binodal we always observe an 
induction time equal to zero, consistent with the nucle- 
ation of small crystallites in the initial stage, and subse- 
quent growth of the crystals. Outside the binodal, the 
induction time is close to zero only at high colloid pack- 
ing fraction. Furthermore, we observe that the fraction 
of crystal is independent of the packing fraction for fixed 
PU, while it decreases for increasing contact energy (3U. 



IV. CONCLUSIONS 

The objective of this work was to systematically 
study the relationship between equilibrium and non- 
equilibrium phases in systems with short-ranged attrac- 
tions. To study the dynamics of the system we carried 
out BD simulations by quenching a homogeneous fluid 
phase to certain state points in the phase diagram. 

First, we find a window of enhanced crystallization 
that is delimited at high contact energy by low diffu- 
sion and at low contact energy (below line A Fig. [5]) by 
a high nucleation barrier. 

Second, we observe that the crystallization is enhanced 
by the metastable gas-liquid phase separation, as demon- 
strated by the fact that the line of spontaneous crystal- 
lization (line A in Fig. [8]) at low density follows closely 
the binodal of the gas-liquid separation (dashed line in 
Fig. [5]). The effect is caused by the two-step process 
crystallization process inside the gas-liquid binodal: the 
system first forms liquid domains, either through nucle- 
ation or spinodal decomposition; subsequently, crystal 
nucleation occurs within the dense liquid domains. This 
process follows the Ostwald rule of stages that states that 
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FIG. 8: (Color online) Sketch of the phase diagram in the 
energy-density representation reporting the location of the 
non-equilibrium phases found in this work. Line A is the line 
of spontaneous crystallization while Line B indicates the tran- 
sition form a cluster phase to a gel. The square insets sketch 
the mesoscopic morphology of the system (clusters, perco- 
lating network (gel) and porous solid from left to right). The 
shading indicate the transition from a more crystalline phase 
(dark red) to a more amorphous one (light yellow). 

the transition from one phase to another may proceed via 
intermediate metastable states. Hence, the effective nu- 
cleation rate of the crystal depends on the nucleation rate 
of the liquid in the gas and the crystal in the liquid phase. 
This two- stage crystallization process was also observed 
in proteins [25|, [26| . 

Last, we find two competing mechanisms of arrest 
of the gas-liquid spinodal decomposition. At low (3U, 
we observe that crystallization starts within the liquid 
branches before the spinodal structure can coarsen into 
macroscopic phase separation. The crystallization slows 
down the coarsening process and arrests the system in 
structures typical of spinodal decomposition. For higher 
(3U, the system becomes amorphous as crystallization is 
inhibited. This result supports a gelation caused by a 
glass transition that meets the liquid binodal and thereby 



arrests the gas- liquid phase separation [22| . We find that 
crystallization and vitrification are both responsible for 
the formation of arrested spinodal patterns (gels), albeit 
at different interaction energies. These two mechanisms 
are also responsible for the formation of two different 
types of long-lived, but transient clusters at low colloid 
density. The clusters are crystalline at low (3U and glassy 
at high PU. At high colloid density, on the right hand 
side of the gas-liquid binodal, the amorphous arrested 
phase is structurally very similar to a (porous) attractive 
glass, while the crystalline arrested phase can be seen as 
a polycrystalline porous material. Fig.[S]shows schemati- 
cally the location of the different type of phases we found 
in relation to the equilibrium phase diagram. 

Finally, we point out that our phase diagram can pro- 
vide a guide to experimentally obtain and control the 
morphology of crystals of colloids or proteins. The crys- 
tallization window can be located experimentally since it 
is delimited at low temperature by an arrested phase with 
liquid-like local structure and at high temperature by a 
mobile fluid phase. Thus, by changing the temperature 
it is possible to control the degree of crystallinity. The 
highest fraction of crystal is obtained for state points 
in the lower interaction strength (higher temperature) 
part of the crystallization window. For increasing inter- 
action strength (decreasing temperature) the crystal frac- 
tion decreases. For the range of packing fractions studied 
in this work {rj c < 0.55), the colloid packing fraction does 
not influence the fraction of crystal or the growth rate, 
but has a strong influence on the morphology of the final 
phase. While at very low rj c clusters form, a branched 
structure is obtained at intermediate packing fractions, 
and porous materials are obtained at high r\ c . 
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